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We present the chemical analysis of 49 giant stars of the globular cluster 
NGC 2419, using medium resolution spectra collected with the multi-object spec- 
trograph DEIMOS@Keck. Previous analysis of this cluster revealed a large dispersion 
^ | in the line strength of the infrared Ca II triplet, suggesting an intrinsic star-to-star 

scatter in its Fe or Ca content. From our analysis, we assess that all the investigated 
stars share the same [Fe/H], [Ca/Fe] and [Ti/Fe] abundance ratios, while a large spread 
in Mg and K abundances is detected. The distribution of [Mg/Fc] is bimodal, with 
~40% of the observed targets having subsolar [Mg/Fe], down to [Mg/Fe]~-1 dex, a 
level of Mg-dcficiency never observed before in globular clusters. It is found that the 
large dispersion in Mg abundances is likely the main origin of the observed dispersion 
of the Ca II triplet lines strengths (that can be erroneously interpreted in terms of 
Fe or Ca abundance scatter) because Mg plays a relevant role in the atmosphere of 
giant stars as an electron donor. A strong depletion in the Mg abundance leads to an 
increase of the line strength of the Ca II triplet, due to the variation in the electronic 
pressure, at a constant Fe and Ca abundance. Finally, we detect an anti-correlation 
between Mg and K abundances, not easily explainable within the framework of the 
current nucleosynthesis models. 

Key words: stars: abundances - stars: atmospheres - stars: evolution - stars: Pop- 
ulation II - (Galaxy:) globular clusters: individual (NGC 2419) 



1 INTRODUCTION 

The old and metal-poor cluster NGC 2419 is by far the 
most luminous globular cluster (GC) residing in the out- 
ermost fringes of the Milky Way (MW) halo, > 10 times 
brighter than any other cluster having Rgc ^ 40 kpc, (see 
iGalleti et all 120071 . for discussion and comparison with re- 
mote clusters in the halo of M31). This unusual feature, 
coupled with a half-light radius much larger than that of 



* Based on data obtained at the W. M. Keck Observatory, which 
is operated as a scientific partnership among the California In- 
stitute of Technology, the University of California, and the Na- 
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ordinary globulars of similar luminosity led several authors 
to the hypothesis that NGC 2419 can be the remnant of 
an originally larger system, like a nucleated dwar f galaxy 
i|Mackev fe van den Berghll2005l ; iBrodie et al.ll201ll ). 

Since the kinematics and the stellar mass function of 
NGC 2419 are fully consistent with a typical GC made 
up of stars and stel lar remnants (see llbata et al.l l2011al : 
iBellazzini et al.ll2012l . and references therein), the key test 
to establish the actual nature of the system is to search 
for inhomogeneities in the abundance of chemical elements 
heavier than Aluminium among member stars. A spread in 
the heavy elements up to the iron-peak group would indicate 
that the system was a site of chemical evolution driven by 
Supernovae (SNe), implying that, at the epoch of SNe ex- 
plosions, the progenitor of NGC 2419 was sufficiently mas- 
sive to retain their highly energetic ejecta (>fewxl0 6 Mq; 
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iBaumgardt. Kroupa fc Parmentierll2008T l , hence it was likely 
a dwarf galaxy. On the other hand, GCs are observed to dis- 
play virtually no spread in heavy elements and large (and 
correlated) spreads in light-elements (Na, O, Mg, Al, in par- 
ticular) , likely associated with a spr ead in He abundance (see 
iGratton. Carretta. B ragaglia 2012, and references therein). 
These signatures are generally believed to trace early chem- 
ical evolution driven by much less energetic polluters than 
SNe, like massive Asymptotic Giant Bran ch stars (AGB 
iD'Antona et all |2002| ; b'Ercole etalll2008l) or Fast Rotat- 
ing Massive Stars CFRMS. see lDecressin et alj|2007al lbl and 
references therein), whose ejecta can be retained in sys- 
tems in the plausible range of mass of proto-GCs. The anti- 
correlation between Na and O abundances is so ubiquitous 
in GCs (and non-existent in the field) that it has been pro- 
posed as a fundamental defining feat ure for globulars (wit h 
respect to open clusters and galaxies ICarretta et alj|2010al ). 

The large d i stance of NGC 2419 (~91 kpc 
|Pi Criscienzo et alj l2011af ) has prevented the detailed 
spectroscopic analysis of a large sample of its stars, that is 
needed to perf orm this chemical te st, until a couple of years 
a g<£J Recentlv lCohen et al.l (|2010l . CIO hereafter) presented 
the analysis of medium-resolution Keck-DEIMOS spectra 
around the Calcium triplet (near 8600 A, CaT) for 43 
Red Giant Branch (RGB) stars of NGC 2419. Using other 
GCs of known [Fe/H] as calibrators and adopting a con- 
stant [Ca/Fe] value they translated the measured pseudo- 
Equivalent Width of CaT and the magnitude difference 
between the stars and the cluster Horizontal Branch (HB, 
V- Vgfl) into [Ca/H] values, as is usually done for [Fe/H] 
fseelRutledge. Hesser fc StetsorJll997l : iBattaglia et~al]|2008l; 
IStarkenburg et ah 2010l . and references therein). CIO found 
a significant spread in the [Ca/H] values derived with this 
method ([Ca/H]caT hereafter): the distribution showed 
a strong peak at [Ca/H]caT — —1.95 and a long tail 
reaching [Ca/H]caT — —1.45. Th is relatively large spread 
in a heavy element (Fe or Ca, see IStarkenburg et alj|20ich 
was interpreted by CIO as providing additional support 
to the hypothesis that NGC 2419 "...is the remnant of a 
dwarf galaxy accreted long ago by the Milky Way...". As a 
consistency check CIO obtained Fe and Ca abundances by 
spectral synthesis of Fe and Ca lines (other than CaT) that 
are present in their spectra. The puzzling result was that, 
while the Ca abundances are in reasonable agreement with 
those from CaT, no significant spread in Iron abundance 
was detected. 

Available photometry is not of much help in settling 
this issue, since in this low metallicity regime the sensitivity 
of optical colours to variations of metallicity is pretty weak. 
Moreover, even if a tiny colour spread is observed along the 
RGB, a spread i n metallicity is not necessar ily the only vi- 
able explanation. iDi Criscienzo et al. I |2011bl ) used exquisite 
HST photometry to show that a small but significant colour 
spread is indeed detected at the base of the RGB, however 
this can be fully accounted for by the Helium spread that 
they invoke to reproduce the complex HB morphology of the 
cluster, within the classical scenario of mult iple populations 
in GCs (|Gratton. Carretta. Br agaglia 2oT3). 
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Figure 1. Metallicity distribution derived from CaT, presented 
as ordinary histograms (lower panel) or as generalised histograms 
(a representation that removes the effects due to the choice of the 
starting point and of the bin width, see lLaird et al. 1988). Note 
that the distribution does not change if only the brightest stars 
are considered (dashed lines). 
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Figure 2. Colour-Magnitude Diagram for the 118 stars included 
in the metallicity distribution from CaT shown in F ig. [T] from 
the photometric catalog assembled by Bcllazzini (2007) (sec Sec- 
tion [2TTJ- Different symbols are adopted for different metallicity 
ranges. A grid of is ochrones o f different metallicity from the BaSTI 
set is overplotted ( Pictrinfcrni et al. 2004), for reference. 



1 IShetrone et al.1 ll200lh provided abundance analysis from high 
resolution spectroscopy for just one star of NGC 2419. 



© 2002 RAS, MNRAS 000, [TV?? 



Chemical composition of NGC2419 3 



Later ICohen. Huang fc Kirbvl (|201ll , Cll hereafter) 
used high-resolution Keck-HIRES spectra to derive the de- 
tailed abundance of several chemical elements in seven 
bright RGB members of NGC 2419. Interestingly, (a) the 
stars studied by Cll do not display any spread in [Fe/H] 
and [Ca/Fe] in excess to what is expected from the ob- 
servational uncertainties, and (b) the abundance pattern of 
these stars is very similar to what is observed in the classi- 
cal metal-poor GC M30 (NGC 7099), residing in the inner 
halo of the MW (D ~ 8 kpc, R G c ^ 7 kpc). One of the 
stars in the Cll sample (SI 131), which is similar to the 
other six in all other aspects, was found to be extremely 
Mg-deficient ([Mg/Fe] = -0.47) and K-rich ([K/Fe]=+1.13), 
a strong anomaly never reported before, at least in Pop-II 
stars (see Cll for a detailed discussion). Hence, at present, 
the results of the chemical test on the nature of NGC 2419 
are still not conclusive. 

Here we report on the chemical analysis of t he medium- 
resolut ion Keck-DEIMOS spectra obtained by llbata et al.l 
|2011al . Ilia hereafter) and used, in that paper, for a 
thorough study o f the kinematics of NGC 2419 (see also 
llbata et alj|2011bh . 



2 OBSERVATIONS 

In this paper we analyse spectra collected with the DEIMOS 
multi-object slit spectrograph mounted on the Keck II 
telescope. This dataset includes spectra of stars located 
along the RGB and observed with the high-resolution 1200 
line/mm grating coupled with the OG550 filter: the spec- 
tra cover the spectral range of 6500-9000 A with a spectral 
resolution of ~1.2 A. Target selection, observations and the 
reduction procedure are described and discussed in full de- 
tail in Ilia. 



2.1 Metallicity from CaT 

[Fe/H]caT was derived by Ilia for the stars observed with 
DEIMOS3 from the pseudo-EW of the three CaT lines (suit- 
ably combined) and V-Vhb (see Ilia for details). Here we 
slightly revise these estimates by (a) adopting the most re- 
cent and reliable value for Vhb , i.e. Vhb = 20.31 ± 0.01 
from Di Criscienzo et al.l (|2011ah instead of Vhb = 20.45 
from lHarrisl ( 1996T ). as done in Ilia, and (b) fully propagat- 
ing all the observational errors into a final error in [Fe/H] CaT 
for each individual star. It is worth noting that there is only 
a tiny difference between the new [Fe/H]caT and those pub- 
lished in Ilia: the new [Fe/H]caT are larger than the older 
ones by 0.037 dex on average, with a r.m.s. of just 0.0004 
dex. 

In Ilia [Fe/H]caT was used as an additional parame- 
ter to clean the kinematic sample from interlopers. While it 
was noted that there was a significant spread in [Fe/H]cuT, 
this result was neither quantitatively assessed nor discussed 
in detail. In Fig. [1] we present the metallicity distribution 



from CaT for the 118 stars from the most conservative 
sample of radial velocity members (their Sample A) hav- 
ing valid [Fe/H]caT in Ilia. The distribution is roughly 
similar to that shown by C10: there is a strong peak at 
[Fe/H]caT — —1.8 and an extended tail (with a possible 
shoulder at [Fe/H]caT — —1.5, not seen in the distribution 
by C10) extending up to [Fe/H]caT — —1.4. Any subtle dif- 
ference in the metallicity zero-point and in the shape of the 
distribution between this work and C10 can be ascribed to 
differences in the adopted calibration and in the sample size 
(118, or 71 stars, vs. 43). Hence we consider the agreement 
of our result with that of C10 as satisfying. 

In the following, to quantitatively estimate the mean 
and the intrinsic spread (a) of a distribution in some ele- 
mental ratio [X/Y] (like, e.g. [Fe/H] or [Ca/Fe]) we adopt 
an algorithm that explores a grid of the ([X/Y],<j) space 
and searches for the couple of parameters that maximises 
the maximum likelihood (ML) function defined as 



ML([X/Y] 
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similar to what is currently done to estimate the mean ve- 
locity and velocity dispersion in kinem atic samples (see, e.g. 
iMartin et aTll2007i : IWalker et ai1l2006l). N is the num ber of 
stars in the sample and a tot = \J& 2 + [X/Y]^ rr . ti , where 
[X/Y]j is the abundance of the individual star and {X/Y] err ,i 
is its associated error. The advantage of this approach is that 
it automatically takes into account the observational errors, 
thus providing estimates of the intrinsic spread with the 
assoc iated uncertainty (derived following IPrvor fc Mevlar] 
1993), a quantity that is especially valuable in the present 
contest. Although the adopted ML function is well suited 
for deriving the parameters of Gaussian distributions, it is 
clearly not a good approximation for the distribution shown 
in Fig. [TJ (but it turns out to be quite appropriate for the 
subsequent applications, see Sect. [3J. Nevertheless, it can 
provide a useful quantitative estimate of the intrinsic spread, 
that is a — 0.17 ± 0.01 dex, similar to the massive globular 
M54 that resi des at the centre of the disrupting Sagittariu s 
dwarf galaxy ilBellazzini et al"1l2008l: [Carretta et aLfeoiOd ). 
A KMM teslPl l| Ashman. Bird fc Zepjl994h . performed both 
on the whole sample or on the sub-sample of the brightest 
stars (V < 19.5), finds that a bimodal model is preferred to 
a unimodal model at the ^ 99.9% confidence level, and the 
relative fraction of stars in the two populations is approxi- 
mately 0.6 and 0.4. 

However, the Colour-Magnitude Diagram (CMD) of the 
target stars shown in Fig. [5] reveals that there are several 
stars that have virtually the same magnitude and color but 
that differ by more than 0.3 dex in [Fe/H]caT- T he adopted 
photo metry is from the catalog assembled by iBellazzim] 
(2007) from various sources; however B, V, I magnitudes 
for most of the target stars (and all the stars we analyse 
with spectral synthesis, below) are take n from the accurate 
absolute photometry by Stetson (2005). To further inves- 
tigate this apparent inconsistency between [Fe/H]caT and 



2 Ilia included in thei r sample also stars from the sample by 
iBaumgardt et al I ^0091 1, for which they took only the radial ve- 
locity estimate obtained by the same authors from HIRES spec- 
tra. 



3 The test compares the likelihood of the data for a single gaus- 
sian model and a double gaussian model and provides the confi- 
de nce level at which the s ingle gaussian model can be rejected, 
see I Ashman. Bird fc Zepl dl994ll . 
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photometry, and having in mind the results by CIO and 
Cll, we decided to perform a full spectral synthesis anal- 
ysis on a subsample of stars of the Ilia sample with high 
Signal-to-Noise Ratio (SNR). This analysis is described and 
discussed in the following sections. 



3 CHEMICAL ANALYSIS 

Among the stars recognised as cluster members by Ilia in 
their Sample A, we selected stars with sufficient SNR (typi- 
cally SNR>40 per pixel) to perform a reliable chemical anal- 
ysis based on spectral synthesis. All the stars have V-band 
magnitude <18.8 and their position in the CMD is shown 
in Fig. [3] Even if the sample is clearly dominated by RGB 
stars, a few Asymptotic Giant Branch (AGB) stars are likely 
included. This is not expected to have any effect on the 
chemical analysis performed below. 



3.1 Atmospheric parameters 

The atmospheric parameters (T e //, surface gravities and 
microturbulent velocities) of the target stars cannot be 
derived directly from the spectra (due to the low spec- 
tral resolution and the small number of available lines) 
and we estimated them from the photometry. Tempera- 
tures are computed by employing the empirical (V — I)p- 
T e ff t ransformation by lAlonso. Arribas fc Martinez-Roger] 
(1999) based on the InfraRed Flux Method. The dered- 
dened color (V — I)o is obtained assumi ng for all the stars a 
color e xcess of E(B-V)= 0.08±0.01 ma g (|Pi Criscienzp et al.l 
2011al) and the extinction law b y iMcCalll (|2004l V The 
Alonso. Arribas fc Martinez-Rogerl l|l999f ) relation for the 
(V — I)o colour is computed in the Johnson photometric 
system, thus, our I Cousin magnitudes have been converted 
into I Johnso n magnitudes by following the prescriptions by 
iBesselll (|l979h . Internal errors in T e f / due to the uncertain- 
ties in photometric data and reddening are of the order of 
about 30-50 K. 

The surface gravities are computed through the Stefan- 
Boltzmann relation, assuming the photometric T e //, and a 
typical mass of 0.75 Mp, according to a n isochrone from the 
BaSTI dataset (|Pietrinferni et al.ll200i ) with age of 13 Gyr, 
Z=0.0003 and alpha-enhanced ([a/Fe]= +0.4 dex) abun- 
dances patterns (see CIO and Cll). For the luminosities, 
we used the dereddened V-band magnitudes with the bolo- 
metric corrections from Al onso, Arribas & Martinez-Roger 
(1999) and the distance modulus of 19.71±0.08 mag 
( Di Criscienzp et al.l l2011aT ). The gravities can be deter- 
mined within ±0.05 dex by propagating the uncertainties 
in T e //, evolutive mass (of the order of 0.05 Mq), distance 
modulus and photometry. 

Finally, the microturbulent velocities are calculated by 
adopting th e rela tion between v t and log g derived by 
iKirbv et al.l |2009r i: this relation provides values of v t of 
~1. 8-2.0 km/s. Formally, the uncertainties in log g trans- 
late into a typical error in v t of 0.01 km/sec. Because the 
majority of the lines analysed here are located in the sat- 
urated/damped portion of the curve of growth, these lines 
are often sensitive to the photospheric velocity fields and 
we adopted a conservative (and more realistic) error of ~0.2 




1 1.5 

V-I 



Figure 3. Color-Magnitude Diagram in the I-(V-I) plane of NGC 
2419 with the spectroscopic targets analysed here mark ed as black 
points . The photometric catalog is that assembled by iBellazzinil 
(2007), in particul ar the source o f the photometry for the spec- 
troscopic targets is IStetsonl ||2005| 1. 



km/s. The adopted atmospheric parameters for all the in- 
vestigated targets are listed in Table 1. 



3.2 Model atmospheres and linelist 

We began the analysis by employing one-dimensional, LTE, 
plane-parallel geometry model atmospheres computed with 
the code ATLAS9 (originally developed by R. L. Kurucz) 
that treats the line opacity through the Opacity Distribution 
Functions (ODFs) method. The ATLAS 9 models are com- 
puted by using the ODFs calculated by ICastelli fc Kuruc3 
(2004) with [M/H]=-2 dex, a-enhanced chemical patterns 
and without the inclusion of the approximate overshooting 
in the calculation of the convective flux. For some stars with 
abundance patterns that differ heavily from those adopted 
in the ODFs computation, specific model atmospheres are 
computed with the code ATLAS 12 that employs the opacity 
sampling method and allows one to calculate mod el atmo- 
spher es with arbitrary chemical composition (see ICastellil 
120051 . for a detailed description of the handling of opacity 
by this code). 

Adopting the photometric atmospheric parameters and 
a guessed metallicity value of [M/H]=-2 dex, we calcu- 
lated synthetic spectra along the entire wavelength range of 
the observed spectra, including all the atomic and molec- 
ular transitions available in the last version of the Ku- 
rucz/Castelli linelist Q. All the synthetic spectra used for 
the analysis have be en generated with th e code SYNTHE 
in its Linux version l|Sbordone et al.ll2004l ; Isbordonell2005h 



4 http:/ /wwwuser. oat. ts.astro.it/castelli/linelists. html 
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and convoluted with a gaussian profile to mimic the spectral 
resolution of the DEIMOS spectra. 

We selected absorption lines that are unblended at the 
DEIMOS resolution and with the atmospheric parameters 
and metallicity of our targets. We included in our linelist 19 
Fe I lines, 5 Ti I lines, the Mg I line at 8806.7 A, the CaT 
line at 8542 and 8662 A and the resonance K I line at 7699 
A (the other K I transition at 7665 A being heavily blended 
with telluric absorption) . Reference Solar values are from 
iGrevesse fe Sauvall jl99sb . 

3.3 Abundance analysis 

The traditional derivation of abundances through the mea- 
sured equivalent widths can be difficult in low resolution 
spectra, due to the line blanketing conditions that make the 
continuum location quite complex. Thus, we resort to the 
comparison with synthetic spectra including all the molec- 
ular and atomic transitions. Each selected line has been 
analysed independently by performing a x 2-mm i m i sa tion 
between the normalised spectrum and a grid of synthetic 
spectra within a spectral window centred around the line 
(in the case of DEIMOS spectra the typical width of the 
fitting window is ~4-5 A). The synthetic spectra were re- 
binned at the pixel-step (0.33 A/pixel) of the observations. 
A preliminary normalisation of all the spectra is performed 
by fitting the entire spectrum with a Chebyshev polynomial. 
For each line, we adjusted locally the normalisation along a 
region of ~60 A around the line. When a first abundance is 
determined from \ 2 minimisation, the spectrum is divided 
by the best-fit synthetic spectrum and the quotient spec- 
tra is fitted with a spline. The resulting spline is used to 
normalise again the observed spectrum and a new x 2 rnin- 
imisati on is performed (the same approach has been adopted 
also bv lShetrone et al.ll2009l ; lKirbv et al.ll20ld) . 

To derive the Ca abundance we employed a slightly 
different methodology, by fitting not the entire line profile 
but only the wings of the two strongest components of the 
CaT, at 8552 and 8662 A. The CaT line at 8498 A has 
been excluded because its wings are marginally sensitive to 
the abundance, at variance to the other two lines that are lo- 
cated in the damped part of the curve of growth and they are 
more sensitive to abundance with respect to saturated lines. 
We excluded the core region of these features because (i) it 
is heavily saturated, (ii) it suffers from relevant departures 
from LTE, and (Hi) its shape depends on the thermal strat- 
ification of the model atmosphere in the outermost atmo- 
spheric layers and basically the canonical one-dimensional 
model atmospheres are not able to describe well the ther- 
mal structure of the upper layers where the line core forms. 

The K resonance line at 7699 A suffers from depar- 
tures from LTE conditions. Unfortunately, large grids of 
corrections as a function of the at mospheric parameters 
are not available for this t ransition. iTakeda et al.l (|2009T ) 
and lAndrievskv et al.l |2010l ) provided abundances of [K/Fe] 
computed by including the departures from LTE for globular 
clusters and extreme metal-poor halo stars. The magnitude 
of their corrections agrees very well and, for the resonance 
li ne used in our work, the correction is quite small (see Fig. 6 
m ITakeda et a"Hl2009h . Because of the small range of metal- 
licity and T e ff of our targets, we decide to apply a unique 
value (-0.30 dex) of the NLTE correction for all the targets. 



Concerning the CaT lines discussed here, NLTE cor- 
rections affiect_c^ily the core of the lines (as pointed 
out by ISpite et al. I 1201 ll . see their Fig. 1, but see also 
IStarkenburg et al* 2010l and references therein), while their 
wings are not affected by the influence of the NLTE and can 
be used to derive reliab l e Ca abundances under the LTE 
assumption. ISpite et al.l (|201ll ) derived Ca abundances of 
metal-poor stars by fitting the wings of the CaT lines with 
both LTE and NLTE profiles, finding that the effect of NLTE 
on the wings of these lines is negligible (less than 0.1 dex). 

Uncertainties in the fitting procedure and in the atmo- 
spheric parameters were taken into account and added in 
quadrature (see Table 1). The uncertainty due to the fitting 
procedure is computed by resorting to MonteCarlo simu- 
lations: for each spectral line, Poissonian noise is injected 
in the best-fit synthetic spectrum (after the re-mapping at 
the pixel-scale of the observed spectrum) and re-analysed 
as described above. One hundred MonteCarlo events are 
computed for each line and the dispersion of the derived 
abundance distribution is taken to be the la uncertainty. 
Typically, the uncertainty in the fitting procedure ranges 
from ±0.09 dex for SNR= 100 to ±0.25 dex for SNR= 40. 
Only for the Ca abundances, the effect of the photon noise 
is quite small (±0.02 dex for SNR= 100 and ±0.05 dex for 
SNR= 40), because of the large number of pixels used in the 
X 2 minimisation. 

For Fe, Ca and Ti (for which we have measures from 
different lines), we computed a weighted average abundance 
for each star by using the error from the MonteCarlo sim- 
ulations as a weight; for the typical internal uncertainty we 
assumed the dispersion of the mean normalised to the root 
mean square of the number of used lines. For K and Mg (for 
which one only line is measured) we take the la level of the 
MonteCarlo distributions as the internal error for each line. 



4 CHEMICAL HOMOGENEITY AND 
INHOMOGENEITY IN NGC 2419 

The derived relative abundances of Fe, Mg, Ti, Ca and K for 
the target stars are reported in Table 1, together with the 
associated uncertainties. The first issue we want to focus on 
is the lack of spread in Fe and Ti abundance in our data. 
Adopting the ML algorithm we find ([Fe/H]) = -2.09±0.01 
with <7 [Pc/H ] = 0.00 ± 0.03, and ([Ti/Fe]) = ±0.29 ± 0.02 
with erpri/Fc] — 0.00 ± 0.04. The agreement with the results 
by Cll is excellent as, with the same algorithm, we obtain 
([Fe/H]) = -2.08 ± 0.05 with <t [Pc/H ] = 0.00 ± 0.08, and 
([Ti/Fe]) = ±0.28 ± 0.05 with cr [Ti/Pc] = 0.00 ± 0.07, from 
their sample of seven starfl 

A small intrinsic spread in [Ca/Fe] is detected, 
([Ca/Fe]) = ±0.46 ± 0.01 with cr [Ca/Fo] = 0.09 ± 0.01, 
while from Cll we get ([Ca/Fe]) = ±0.14 ± 0.06 with 
""[Ca/Fc] = 0.00 ±0.07. The non-null spread in [Ca/Fe] found 
in our data may be partially due to an underestimate of the 
errors on the Ca abundance we derive from the wings of 
the CaT lines, since these lines are in a different part of the 



5 We have considered here the abundances derived by Cll from 
Till lines. From Til ([Ti/Fe]) = +0.12 ± 0.02 with <T[xi/Fe] = 
0.00 ±0.07 
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curve of growth with respect to all the other lines consid- 
ered here. Correlation between uncertainties in Fe and Ca 
abundances may also play a role, as [Ca/Fe] appears to cor- 
relate with [Fe/H] even if the spread in [Fe/H] is not statis- 
tically significant, while the correlation between [Ca/H] and 
[Fe/H] is much weaker, if any. This view is supported by 
the fact that we do not detect an intrinsic spread in [Ca/H]: 
([Ca/H]) = -1.63 ± 0.01 with <r [G a/H] = 0.00 ± 0.02. We 
conclude that the spread in Calcium abundance is absent or 
very small. 

Contrary to Fe, Ca and Ti, large star-to-star spreads 
in the Mg and K content are detected among the stars of 
NGC 2419: ([Mg/Fe]) = +0.05±0.08 with cr [Mg/Fc] = 0.56± 
0.06 and ([K/Fe]) = +0.92±0.08 with <t [k /f c ] = +0.51±0.06. 
As we shall see in detail below, a single gaussian model 
in not appropriate for the observed [Mg/Fe] and [K/Fe] 
distribution, since these are clearly bimodal. The average 
and intrinsic spread for [Mg/Fe] are in reasonable agree- 
ment with the Cll results ({[Mg/Fe]) = +0.29 ± 0.13 with 
® [Mg/Fe] = 0.30 ± 0.10; see Sect.O below). 

Cll provides K abundance only for the anomalous star 
SI 131. The discussion of the K abundances is deferred to 
Sect. U 

The star Iball_83 was also observed by Cll (star S1209 
in their sample): our T e ft is lower than that derived by Cll 
by 200 K and this difference arises mainly from the differ- 
ent colour- T e ff transformations and colour excess employed 
in the two studies. We found that the differences with re- 
spect to Cll are basically ascribable to the difference in T e ff; 
only for Ca is the discrepancy large. Taking into account the 
corresponding iron abundances, our Ca abundance is higher 
than that found by Cll by ~0.3 dex, and this difference 
is only partially explainable with the employed T e g scales. 
We can suppose that the residual difference is due to the 
adopted Ca lines. 

4.1 Metallicity from CaT in Mg-deficient 
atmosphere 

The large intrinsic spreads in [Ca/H]caT and in [Fe/H]caT 
reported by C10 and Ilia, respectively, and found again in 
Sect. 12. ll above, requires an additional discussion, since our 
spectral synthesis analysis, in fact, confirmed the results by 
Cll, i.e. there is no intrinsic scatter in Iron abundance and 
just a small (if any) spread in Calcium. 

The observed large variance of the CaT strength is not 
due to an Fe or Ca abundance spread but might be explained 
in light of the large range covered by the Mg abundance. In 
fact, Mg plays a relevant role in the opacity of cool stars, be- 
ing one of the most important electron donors, contributing 
to the formation of the H~ ions, the main source of opacity 
in these stars together with Fe, Al and the other a-elements 
(the relative contribution of each element depends on the 
temperature, gravity and optical depth). A depletion of Mg 
implies a lower number of free electrons available to form H~ 
ions, leading to a decrease of the H~ opacity. The strength 
of the CaT lines in creases, decreasing th e electron pressure 
(as pointed out bv lShetrone et all 2009). Thus, a large de- 
pletion of the Mg abundance produces an increase of the 
strength of the CaT lines at fixed Ca and Fe abundance. 

Fig.|4]shows the comparison of synthetic spectra around 
the Mg I line (first panel) and the three CaT lines computed 
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Figure 4. Comparison between the synthetic spectra around the 
Mg I lines at 8806.7 A (upper-left panel) and the three CaT 
lines computed by employing ATLAS 12 model atmospheres with 
T e// =4200 K, log g=0.6, [M/H]=-2 dex, ^-enhanced patterns, 
but with [Mg/Fe]=+0.4 dex (black spectrum) and [Mg/Fe] =-1.0 
dex (red spectrum). The synthetic spectra are at the same spec- 
tral resolution as the DEIMOS spectra. 



by adopting ATLAS 12 model atmospheres with T e ff = 4200 
K, log g= 0.6 and [M/H]=-2 dex, and with [Mg/Fe]=+0.4 
dex and [Mg/Fe] =-1.0 dex (black and red line respectively). 
In the latter case, the strength of the CaT lines significantly 
increases with respect to the synthetic with [Mg/Fe] =+0.4 
dex (note that the lines of elements other than Mg and Ca, 
that are also visible in the spectra, are identical in the two 
cases). We stress that the two synthetic spectra are iden- 
tical in terms of overall metallicity, individual abundance 
patterns (in particular the same Fe and Ca abundances) 
and atmospheric parameter, and the only difference is in 
the Mg abundance. Fig. [S] illustrates the behaviour of the 
electronic pressure and of the temperature as a function of 
the Rosseland opacity for the two ATLAS 12 model atmo- 
spheres, where the differences in terms of T e ff and P e are 
appreciable (in particular in the outermost layers where the 
core of the CaT lines forms and in the deepest regions where 
the wings of CaT form). 

As a confirmation of this effect, we detect a clear 
anti-correlation between our [Mg/Fe] abundance ratio and 
[Fe/H] CaT (Fig. The iron-rich stars (with [Fe/H] CaT >- 
1.7) dex have [Mg/Fe] <0 dex, and the bimodality in 
[Fe/H]caT noted in Sect. [2~T1 reflects the bimodal distribu- 
tion of Mg abundance (see below) . 

The link between the Mg abundance and the intensity 
of the CaT lines is appreciable also in Fig. [7] that compares 
portions of the spectra for the two target stars Iball_28 and 
Iball_3 (red and black line respectively), that have very 
similar T e ff and log g (see Table 1) and different [Fe/H]caT 
(by more than 0.4 dex, but indistinguishable Iron abundance 
as derived from Fe lines). The difference in the line strength 
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Figure 5. Comparison between the two ATLAS12 models de- 
scribed in the caption of Fig. [4] Upper panel: electronic pressure 
as a function of the Rosseland optical depth. Lower panel: tem- 
perature as a function of the Rosseland optical depth. The grey 
regions indicate the regions of line formation for the core and the 
wings of the Ca II triplet 
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Figure 7. Comparison between the spectra of the stars Iball_28 
(red line) and Iball_3 (black line) around the CaT line at 8662 
A (upper panel) and around the Mg I line at 8806.7 A. Note 
that the two spectra have very similar atmospheric parameters 
but different [Fe/H] based on the CaT lines. In particular they 
have T ef f = 4201 °K and T eff = 4224 °K, log g=0.60 and 
log g=0.62, respectively, and the same value of microturbulent 
velocity, vt = 1.99 km s _1 . 
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Figure 6. Behaviour of the [Mg/Fe] abundance ratio as a func- 
tion of the metallicity derived from the CaT lines. The [Mg/Fe] 
abundances do not include NLTE corrections: as discussed in Sec- 
tion [5] the NLTE corrections are always negative (thus lowering 
the [Mg/Fe] ratios) and their magnitudes are a function of the 
Mg abundance: the most Mg-rich stars can be overestimated by 
0.2-0.3 dex, while the most Mg-poor stars are basically unaffected 
by departures from NLTE. 



of the CaT line at 8662 A is evident (upper panel), as well 
as the difference in the Mg I line in the lower panel. On 
the other hand, the Fe I lines visible in the spectra are very 
similar to each other. 

Note that the Ca abundances that we derived from spec- 
tral synthesis of the wings of the CaT lines take into account 
the effect of the Mg abundance: the wings were fitted by us- 
ing synthetic spectra computed with the corresponding Mg 
abundance derived from our analysis 0. When we repeated 
the fitting procedure by using only synthetic spectra with 
[Mg/Fe] =+0.4 dex (the value of the employed ODFs), the 
Ca abundance of the Mg-deficient stars are over-estimated 
by about 0.2-0.3 dex. This difference is smaller than the Ca 
abundance range found by C10 (spanning about 0.6-0.7 dex) 
because the measurement of the equivalent widths of these 
lines includes also the contribution of the line core, that is 
highly affected by the variation of the electronic density (see 
Fig. [4]) but less sensitive to the Ca abundance. 

The Mg depletion observed among the stars of the clus- 
ter might explain both the lar ge distributions of [C a/H] CaT 
by C10 and of [Fe/H] CaT by llbata et all i|2011al ). As ad- 
ditional confirmation of this interpretation, the Mg-poor 
star S1131 discussed by Cll has also one of the highest 
[Ca/H]caT values in the whole C10 sample. 



6 For sake of completeness we repeated the analysis of all the 
stars with [Mg/Fe]<-0.2 dex by using ATLAS 12 model atmo- 
spheres computed with the corresponding abundances. Basically, 
the difference with respect to the analysis performed with the 
ATLAS 9 models is of a few hundredths of dex. 
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Finally, we stress that the strength of CaT lines depends 
on the global budget of the free electrons, thus, in addition 
to Mg, other elements that are electron donors in the at- 
mospheres of cool stars can contribute to shape the CaT 
lines. In particular, the Al abundance deserves a note of 
caution, since in ordinary GCs [Al/Fe] anti- correlates with 
[Mg/Fe] (iGratton, Carretta, Bragaglial I2012T ). In principle, 
an enhancement of Al abundance leads to an increase of 
the free electrons, that could counterbalance those lacking 
because of Mg depletion. This may be a source of concern 
for the interpretation of the abundance spreads derived from 
CaT presented above, since, unfortunately, we lack Al abun- 
dance for our target stars. We investigated the effect of dif- 
ferent Al abundances (in the range 0.0 ^[Al/Fe]^ +1.5) 
coupled with a Mg depletion ([Mg/Fe]=-1) on the strength 
of CaT lines, using suitable ATLAS12 model atmospheres. 
We found that the strength of CaT decreases appreciably 
only for [Al/Fe] ^1 dex, thus mitigating the effect of the 
Mg depletion only in this extreme regime. It is interesting 
to note that Cll finds an Al abundance much lower than 
this for the only Mg-deficient star in their sample, S1131 
( [Al/Fe] — +0.45). Assuming that this value is representative 
for the Al abundance in Mg-deficient stars of NGC 2419, the 
effect on the CaT lines should be negligible. 



5 MAGNESIUM DEPLETION 

The stars in our sample appear to belong to two groups, 
according to their Mg abundance (see Fig. [5] and Fig. [5]). 
A first group is more or less in the range c overed by stars 
in ot her GCs, 0.0 < [Mg/Fe] < +0.8 dex |Carretta et all 
2009). Consequently we will dub these stars Mg-normal. The 
second group is made of Mg-deficient stars having — 1.4 < 
[Mg/Fe] < 0.0 dex, a range usually not reached by stars in 
Galactic globulars (see below). The KMM test confirms that 
the [Mg/Fe] distribution is bimodal with a confidence level 
> 99.9%. A fraction of 62% of the stars are attributed to the 
Mg-normal population, having ([Mg/Fe]) = +0.47 dex, and 
the remaining 38% to the Mg-deficient population, having 
([Mg/Fe]) = -0.62 dex. 

Intrinsic dispersions in the Mg abundance are usually 
observed in GCs and explained in the framework of self- 
enrichment processes, occurring in the early stages of GC 
evolution. The majority of the Milky Way GC stars for 
which Mg abundances are available covers a range from +0.2 
to +0.5 dex, the only known cluster having Mg-deficient 
stars is NGC 2808, where three RGB stars w ith [Mg/Fe] be- 
tween -0.3 and -0.1 dex have been detected l|Carretta et al.l 
I2009T ). Note that NGC 2808 is characterised by a huge 
Na-O anti-correlation, with the presence of super O-poor 
stars, and (likely ) a wide and multi-modal distribution in 
He abundances (iD'Antona et all 120051 ; iPiotto et all 120071 ; 
iBragaglia et alJl201Ch . Also, two RGB stars with [Mg/Fe] =- 
0.21 and -0.31 dex have been meas ured in the Large Mag - 
ellanic Cloud old GC NGC 1786 (|Mucciarelli et al l I2009T I 
and both the stars share a very large (<-0.4 dex) deple- 
tion of [O/Fe]. Hence Mg-deficient stars are rare in GCs 
and seem to be linked to the hypothesised second genera- 
tion of stars formed from the e.je cta of massive AGB st ars 
|D'Ercole et all 120081 ) or FRMS l|Decressin et all l2007bl ) of 
the first generation, and characterised by high helium and 



sodium content and strong deple tion of O and Mg (see 
IGratton. Carretta. Bragaglial [20121 and references therein). 

Recent an alyses of the bimodal HB morphology 
of N GC 2419 (|Dalessandro et all 120081 ; ISandauist fc Hessl 
2008) show that a relevant fraction (~30%) of the HB pop- 
ulation is composed of extremely hot stars which have been 
tentatively interpreted as seco nd-generation stars with hig h 
helium content (Y= 0.42) by |Pi Criscienzo et all (|2011bl ). 
The fraction of these extreme HB/Blue Hook stars is fairly 
similar to the fraction of Mg-deficient stars observed in our 
sample, and the postulated bimodality in the He content 
is intriguingly reminiscent of the bimodal distribution of 
[Mg/Fe] detected here. 

However the range of [Mg/Fe] covered by our stars is 
much larger than anything ever seen before in any stellar 
system. Hence, while we will come back to discuss the pos- 
sible physical origin of this feature in Sect. [S] it is worth 
discussing briefly some details of the analysis that may lead 
to erroneous estimates of [Mg/Fe]. 

In Fig. [7] it is clearly shown that large differences in 
the strength of Mg lines are observed in the spectra of stars 
with the same atmospheric parameter^!, hence some factor 
other than a variation in these parameters should be respon- 
sible of the large star-to-star differences in Mg lines. This is 
confirmed also by the upper panels of Fig. [8] showing that 
there is no discernible trend between Mg abundance and 
T e ff and log g. Systematics in the estimates of the micro- 
turbulent velocity (the only atmospheric parameter that we 
cannot derive directly from photometry) does not seem able 
to produce spurious underestimates of the Mg abundance. 
In fact, the Mg lines for stars of similar T e ff and log g are 
distributed in both the linear and saturated portion of the 
curve of growth: the weakest lines (those with anomalous low 
Mg abundances) are those basically unaffected by the veloc- 
ity fields. Changes in microturbulent velocities will change 
the Mg abundance for the Mg-rich stars but will have a small 
or negligible impact on the Mg-poor stars. 

Furthermore, the lack of trend of [Mg/Fe] ratio as a 
function of T e f j and log g suggests that NLTE effects play a 
negligible role and they are n ot responsib l e for t he observed 
Mg range. As pointed out bv lMerle et~"ai1 f|201lh . the NLTE 
corrections (in the sense Anlte-Alte) of the Mg line for 
a giant star with [M/H]=-2 dex and a-enhanced chemical 
mixture are negative, hence the application of the correction 
would lead to even lower Mg abundance for Mg-deficient 
stars. 

In order to check the magnitude of the NLTE correc- 
tions for this Mg line, we calculated the departures from 
LTE by using a mod ified version of the NLTE radiative 
transfer code MULTI (ICarlssonl ll986) and a -poor MARCS 
model atmospheres (|Gustafsson et al.ll2008l ) with the typi- 
cal parameters of our targets: T e g= 4250 K, log<7= 1.5 and 
[Fe/H]=-2. We used a simplified model atom of Mg I with 
15 energy levels, 26 radiative bound-bound transitions and 
91 electron collisional bo und-bound transi tions (the atomic 
data are the ones used in lMerle et al.ll2~01ll ). We assumed no 
contribution from inelastic collision with hydrogen: this will 



7 The considered stars differ by just 23 °K in the adopted T e ff, 
by 0.02 dex in log g, and by 0.0 km s — 1 in vt, that is significantly 
less than the uncertainties in each parameter. 
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produce upper limits in NLTE abundance corrections. The 
magnitude of the corrections decreases for lower Mg abun- 
dances, because the line is formed progressively deeper in the 
photosphere, reducing the contribution of the core (where 
the deviations from LTE are more relevant): for [Mg/Fe] =- 
1 dex the NLTE correction is <0.1 dex and the line profile 
is basically the same under both LTE and NLTE assump- 
tions, while the maximum correction (-0.35 dex) is for high 
Mg abundances. Thus, we can exclude that the very low 
[Mg/Fe] values are due to NLTE effects. On the other hand, 
proper NLTE corrections may lower the Mg abundance of 
the most Mg-rich stars in our sample, thus leading to an 
upper limit in the [Mg/Fe] distribution of NGC 2419 more 
similar to that observed in other Galactic GCs of compa- 
rable metallicity (see Fig. II Ut below). Lacking a sufficiently 
fine grid of NLTE corrections we are forced to neglect this 
factor: the reader must be aware that the Mg abundance 
of stars having [Mg/Fe] > +0.5 can be overestimated by 
0.2-0.3 dex. 

As a sanity check of the Mg indicator used here, we 
analysed the DEIMOS spectrum of the star S1131 (kindly 
provided by E. Kirby and J. Cohen) that Cll identify as 
Mg-poor through the measurements of other Mg lines. When 
the atmospheric parameters by Cll are adopted, we derive 
[Mg/Fe]=-0.35±0.13, in perfect agreement with the HIRES 
results. We note that our T e // scale is slightly hotter than 
that used by Cll (due to the adopted reddening value and 
T e //-colour transformations). Still, if we derive the atmo- 
spheric parameters for S1131 following the same procedure 
adopted for the stars in our sample, as described in Section 
3.1, we obtain [Mg/Fe]=-0. 27+0.13, again in good agree- 
ment with Cll. We stress that the [Mg/Fe] estimates by 
Cll and by us for S1131 are based on different Mg lines and 
the slightly higher abundance derived from the Mg line at 
8806.7 A can be easily explained by taking into account the 
NLTE effects (note however that the two values are fully 
consistent within the errors, the actual difference being 0.08 
dex). 

Finally, we checked the impact of a strong helium en- 
hancement and compare synthetic spectra computed by us- 
ing ATLAS 12 model atmospheres with normal Y abundance 
(Y=0.24) and with Y= +0.4. The net effect of a high Y 
value is to make the lines deeper (especially for the strong 
features) but this effect is very small (typically less than 
0.5% of the line depth computed with normal Y) and to- 
tally negligible in the analysis of real spectra. Note that the 
assumption of Y= +0.4 does not modify the line profile of 
the Mg I line at 8806.7 A, thus we rule out that the strong 
depletion of Mg abundance is due to the effect of a very 
strong He enhancement. 

Unfortunately, we cannot measure other elements in- 
volved in the chemic al anti-correlations typical of GCs 
l|Carretta et alJfeoiOal ). O and Al cannot be detected at this 
spectral resolution, and the only Na transitions available 
in our spectral range are the doublet at 8183-94 A, that 
is heavily contaminated by telluric lines (we highlight that 
the cleaning of telluric lines at this resolution can be quite 
imprecise). We cannot assess if the Mg depletion is linked 
to variations of the other light elements. However, the very 
large Mg distribution of NGC 2419 could be the clue of a 
peculiar star formation history of the cluster and indicate 
extremely extended anti-correlation patterns. 



• 

•• • + 
• • • • • 
• 

•• 


• • *• • - 
• 

• • • • • 

— • — 
• • 

i t i i i t i i r t i t 


• 


i i i i i i i i i i i i 

• V 

• r • 

• • -A 

z • 



4000 4500 0.5 1 



T„ H (K) logg 

Figure 8. Behaviour of [Mg/Fc] (upper panels) and [K/Fe] (lower 
panels) abundance ratios as a function of T e f e (left panels) and 
log g (right panels). 
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Figure 9. Behaviour of [K/Fe] abundance ratios as a function 
of [Mg/Fe] ratios. Black triangles are the stars analysed by Cll. 
The histograms of the [Mg/Fe] and [K/Fe] distributions are also 
plotted. 

6 A MG-K ANTICORRELATION? 

Fig. [5] displays the clear and significant anti-correlation be- 
tween [Mg/Fe] and [K/Fe] that occurs in our sample. Cll 
dedicated a detailed description to their star S1131, labelled 
as "a star with a peculiar abundance pattern" , characterised 
by Mg depletion ([Mg/Fe] =-0.47 dex) and K enhancement 
([K/Fe]= +1.13 dex). These abundances nicely match those 



© 2002 RAS, MNRAS 000, [L>?? 



10 Mucciarelli et al. 



obtained in our analysis (the stars analysed by Cll are plot- 
ted for reference in Fig. [9] as black triangles). Taking advan- 
tage of our sample of 49 stars, we can assess that this star 
is not an isolated case, but it belongs to a sub-population of 
stars in NGC 2419 characterised by sub-solar [Mg/Fe] ratios 
and high (>1 dex) [K/Fe] ratios. 

We investigated some technical reasons that may be 
able to mimic the observed anti-correlation. 

(i) Contamination by telluric lines. The used K line is lo- 
cated near the red side of the telluric absorption A band, 
hence it may be prone to some contamination. We checked 
accurately for possible blending with telluric lines, by us- 
ing as template the DEIMOS spectrum of a rapidly-rotating 
white dwarf. In all the spectra the K line is free from telluric 
lines. The Mg line is far from any telluric feature. 

(ii) Errors in the atmospheric parameters. Both the 
K and Mg abundance are virtually insensitive to the 
adopted gravity (at a level of ±0.01 dex for a variation of 
Slog g=±0.1). The abundances are more sensitive to the 
T e ff. a variation of ±50 K translates into a variation of 
the Mg number density of ±0.06 dex and of the K number 
density of ±0.1 dex. The errors in T e // can lead a correla- 
tion between the Mg and K abundances but are not able 
to mimic a spurious anti-correlation. In our analysis v t is 
not constrained spectroscopically (as is usually done) but 
we adopted an empirical relation between vt and log g. A 
variation of ±0.2 km/s leads to a variation of =p0.08 dex and 
^0.14 dex in Mg and K abundance respectively. Also in this 
case the errors run in the same direction and are not able 
to introduce an anti-correlation; 

(iii) Temperature scale. We checked the impact of a differ- 
ent T e ff scale, assuming the temperatures inferred by pro- 
jecting the position of each target in the CMD on the best-fit 
isochrone (see Section 3.1). The average difference between 
the two scales is T ^° s chrone -T^] naom =-%7 K (<r= 52 K). 
The chemical analysis with this T e // scale does not signif- 
icantly change the observed dispersions both in Mg and K 
abundances, and does not erase the Mg-K anti-correlation. 

(iv) Departures from LTE. We applied a unique NLTE 
correction to all the targets, following the work by 
iTakeda et al] |2009l ). The NLTE corrections are usually a 
function of the metallicity and the atmospheric parameters: 
we note that stars in our sample with the same parameters 
exhibit different K abundances, suggesting that our NLTE 
correction is not the cause of the large dispersion in the 
abundance of Potassium. 

The two elements are mainly produced by massive stars: 
K comes from hydrostatic oxygen shell burning and explo- 
sive oxygen burning, while Mg is produced by hydrostatic 
carbon burning and explosive neon burning. However, it is 
noteworthy that the current yields for K are inadequate to 
repr oduce the observed c hemical patterns in the Milky Way 
fsee lRomano et al"1l2010T l. Also, Mg is destroyed through the 
Mg-Al cycle of the hot CNO chain in AGB stars. However, 
production of K through the proton-capture on Argon nuclei 
could occur in the nucleosynthesis chains of massive AGB 
stars. Unfortunately the uncertainties related to the cross- 
sections of these reactions prevent any reliable predictions; 
preliminary computations indicate that a simultaneous vari- 
ation of the initial Argon abundance and of the cross-section 
of the proton- Argon capture can produce ejecta as enriched 



in K as we observe in the atmosphere of the Mg-deficient 
star s (P.Ventura and F. D'Antona, private communication). 

ITakeda et al.l (|2009h identify two K-rich stars in GCs 
M 4 and M 13, coupled with emission components along the 
H a line profile. These findings suggest an increase of the 
velocity fields in the upper photospheric layers where the 
core of the K I line at 7699 A forms. On the other hand, Cll 
ruled out this hypothesis, because some stars with normal K 
abundances exhibit H a emissions larger than those observed 
in the K-rich star S1131. Unfortunately, the setup of our 
spectra does not include the H a line and we cannot check 
this possibility. In any case, a rise of the turbulent velocity 
field in the outermost layers of the photosphere should affect 
also the core of the Mg I at 8806.7 A in a similar fashion. 

Finally, we note that the two employed transitions have 
very different excitational potentials (0 and 4.34 eV for K 
and Mg line respectively), thus they formed in different 
photospheric layers. In fact, the formation of low-x lines is 
favoured in external layers, due to the decrease of the tem- 
perature with depth, while high-x transitions occur prefer- 
entially in the deepest photospheric regions. In particular, 
the wings of the two lines form at similar depths, while the 
formation depth of the core is shifted towa rd low r, th us de- 
creasing x (f° r comparison see Fig. 13.4 in lGravl [l992'). The 
K line forms in the outermost photospheric region where 
the thermal fluctuations can be relevant (in contrast to the 
Mg line that forms deeply). Also, temporal variations of the 
thermal structure of the photosphere cannot be excluded 
and they could introduce an additional scatter in K abun- 
dances but they do not totally explain the anti-correlation 
with Mg. In any case, the observed spread in Potassium re- 
quires an independent confirmation using other lines, as well 
as a search for the Mg-K anti-correlation in other clusters. 



7 A COMMENT ABOUT CAT 

The CaT represents a widely used indicator of metallic- 
ity, especially for distant stellar populations, for which high 
SNR, high resolution spectra are not easy to achieve. Basi- 
cally, the combined EW of the CaT lines is a function of the 
met allicity and it is interpret ed in terms of [Ca/H] or [Fe/H] 
(see lStarkenburg et al.ll2010l . and references therein). 

From a theoretical point of view the strength of the 
CaT lines is sensitive to the abundance of Ca (of course) 
and Fe, but also to the abundance of those elements that af- 
fect the H~ continuum opacity, through their contribution 
to the electronic density. Basically, the calibrations of the 
CaT EWs are referred only to cases where all a-elements 
abundance patterns are described by an unique value: for 
instance, stars with the same level of enhancement of each 
a-element. The [a/Fe] value sets not only the abundance 
of each a-element but also the free electron density, which 
has a relevant impact in the continuum opacity of cool stars 
(as discussed above). Spurious cases, where the a-elements 
have largely different patterns to each other, are not taken 
into account in the usual CaT calibration, but, as demon- 
strated by the case of NGC 2419, they may lead to a wrong 
interpretation of the CaT strengths. In fact, the CaT can 
be used as a diagnostic of the metallicity only for stars with 
canonical chemical composition, while the method fails with 
stars with a peculiar or exotic chemical mixture, because the 
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effective global budget of the free electrons is not taken into 
account. 

Mg deficiencies as extreme as those observed here seem I ' r 

to be very rare in nearby stellar systems (see Fig. 1101 be- 
low, and discussion in Cll, Sect.[5]and Sect. [8]). Therefore it 
is unlikely that their metallicity distributions derived from 1 
CaT are seriously biased by this effect. However, in princi- 
ple, observed dispersions of CaT in GCs (or dwarf galaxies) . '-^ii^i^iTf Af* 
can be due to (i) an intrinsic variation of Fe (under the as- 
sumption that all the stars share the same a content), (ii) £ 
an intrinsic variation of Ca (under the assumption that the ~m 
abundances of the electron donors do not vary from star ■— ' 
to star), but also to (Hi) an intrinsic dispersion of one (or 
more) of those elements that are electron donors. Hence, in 
deriving metallicities of individual stars from CaT lines, it 

seems advisable to counter-check also the abundance of Mg ~ 1 ° 

(which should be feasible to undertake in most cases, since o Q 
the Mg I line at 8806.7 A lies very close to the CaT), to 

validate the adopted calibration and avoid possible spurious _ 2 
effects as those described in Sect. 14.11 [Fe/H] 




8 SUMMARY 

The main results derived from the analysis of 49 giant stars 
in NGC 2419 are summarised as follows: 

• All the stars share the same iron content, with an av- 
erage iron abundance of [Fe/H]=-2.09±0.02 dex (cr= 0.11 
dex), where the observed dispersion is fully compatible with 
the uncertainties. Also, [Ca/Fe] and [Ti/Fe] turn out to be 
homogeneous; 

• NGC 2419 exhibits a large dispersion in the Mg abun- 
dance, reaching values of [Mg/Fe]~-1 dex (unusual for GC 
stars). The large spread of Mg is likely the main origin of 
the observed dispersion of the CaT lines strength, previously 
interpreted as intrinsic dispersion in Ca or Fe. In fact, a Mg 
depletion leads to an increase of the equivalent widths of the 
CaT lines (at a constant Ca abundance). This effect is con- 
firmed by the fact that the iron content inferred from the Ca 
triplet lines clearly anti-correlates with the Mg abundances. 
However, we bear in mind that the strength of the CaT lines 
suffers from the abundances of all the elements that are elec- 
tron donors and the correct line profile for these lines should 
be derived by taking into account the global budget of the 
free electrons; 

• The [Mg/Fe] distribution is bimodal, with about 40% 
of the stars having sub-solar [Mg/Fe] ab undance ratio. This 
fraction is similar to that suggested bv lDi Criscienzo et al.l 
(|2011bh for the extreme population stars with initial helium 
abundance of Y~+0.4, according to the Horizontal Branch 
morphology of the cluster; 

• a very large spread in K content is detected among 
the stars of NGC 2419, spanning from solar values up 
to [K/Fe]~2 dex, with a bimodal distribution. A clearcut 
anti-correlation between [Mg/Fe] and [K/Fe] is observed, in 
agreement with the results already found by Cll for the star 
S1131, that shows an unusual depletion of Mg coupled with 
a strong enhancement of K. 

Fig. [10] provides a direct illustration of the extremely 
unusual abundance pattern of NGC 2419. The spread in 



Figure 10. [Mg/Fc] vs. [Fc/H] for a) open circles: stars of 
NGC 2419 from the present study, b ) triangles: stars of NGC 2419 
from Cll, c) s mall dots: stars in the dwarf spheroidal satellites of 
the MW from lKirby et al.l j201Ch d) X symbols: stars of various 
Galactic GCs from lCarretta et al.l (I2009T) . Stars from NGC 7078 
and NGC 2808 are plotted in blue and red, respectively, for refer- 
ence, since they displ ay the most Mg-deficient stars in the whole 
ICarretta et al.l i2009T ) sample. 

[Mg/Fe] is unrivalled both in GCs and in dwarf galaxies. A 
strong depletion of Mg is generally interpreted as due to a 
significant contribution of ejecta from S NIa to the gas from 
which Mg-poor stars are formed (see iTsuiimoto fe Bekkil 
|2012| . and references therein). Indeed, Fig. llOl shows that Mg- 
depletion with increasing metallicity is a common feature 
in dwarf spheroidal (dSph) galaxie s and that [Mg/Fe] < 0.0 
stars are not rare in those systems. iLemasle et alj (|2012t ) re- 
ports that the two most metal-rich stars in their sample of 
red giants in the Carina dSph have —0.6 <[Mg/Fe]< —1.0, 
not far from the most Mg-deficient stars in our NGC 2419 
sample. However, the Mg depletion in dSph (as in any other 
galaxies studied until now) is always coupled with an in- 
crease of [Fe/Hfl since the whole effect is due to SNIa en- 
riching the interstellar medium with material that is Fe- 
rich and poor in a elements, thus reducing the [Mg/Fe] ra- 
tio with respect to the pattern previously set up by SNII 
(see iTolstov. Hill, fe Tosil 120091 . for a recent review). This 
clearly does not occur in NGC 2419 where stars in the 
range —1.4 <[Mg/Fe]< +0.8 are indistinguishable in terms 
of [Fe/H]. 

This fundamental observational fact, coupled with the 
analogy of the bi-modalities in [Mg/Fe] and HB morphol- 
ogy recalled above, suggest as more likely the possibility 
that the observed pattern was produced by the same pro- 
cesses that cause the Na-O anti-correlations, and other sig- 

8 As well as a star forma tion history lasting for a few Gyrs 
I Tolstoy, Hill, fe Tosil |2009|) , clearly not observed in NGC 2419 
I Pi Criscienzo et a.l.ll2011bl ; iBellazzini et alj|2012l) . 
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natures of the early chemical enrichment that are peculiar 
to GC s |Carretta et al.ll2010al ; lGratton. Carretta. Bragaglial 
|2012| ). Further support to this hypothesis is provided by the 
fact that Mg-deficient stars lie systematically to the red of 
Mg-normal stars along the RGB, in the V, U-V Colour- 
Magnitude Diagram (Lardo et al., in preparation). In the 
framework of anti-correlations in GCs, Mg-deficient stars 
should correspond to Na-rich (and O-poor) stars and it is 
a generally observed characteristic of GCs to have Na-rich 
RGB stars that are redder than their Na -poor counterparts 
in CMDs that include the U-band (see lLardo et al.l l201ll . 
for discussion and references). The colour difference is due 
to variations in the strength of CN and CH features at 
wavelengths lower than ~ 4000 A, that are driven by vari- 
ation of C and N abundance, that, in turn, are correlated 
with the abundance of other light elements tracing the self- 
enric hment process in GCs (like Na and O. ISbordone et all 

Eon]). 

Even if the observed Mg spread can be attributed to 
this kind of processes, the reason for the extreme behaviour 
of this cluster, as well as the origin of the Mg-K anti- 
correlation, remain unclear. This latter feature, if confirmed, 
may shed new light into the whole process of GC formation, 
possibly providing a tool that may help to d iscriminate be- 
tween competing models of s elf-enrichment l| D'Ercole et al.l 
120081 : iDecressin etafl l2007al lbl). An obvious singular char- 
acteristic of NGC 2419 is its distance from the centre of 
the Galaxy. Is it possible that a low degree of interaction 
with - for instance - the gas rich disc of the Milky Way has 
favoured the retention of enriched material in this cluster. 
iBellazzini et al.l (|2012h found some indirect evidence for the 
cluster having a significant larger t o tal ma ss in the past, 
as postulated bv lDi Criscienzo et al.l (|2011bh on completely 
a different basis: also a large total stellar mass may have 
played a role (but this is a general feature of models of GC 
formation accounting for the presence of multiple popula- 
tions). In any case, this mysterious cluster does not cease 
to reveal new interesting features that makes it especially 
worthy of further study. 
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Table 1. Identification numbers (from Ilia), temperatures, gravities, microturbulent velocities and abundance ratios for 
the observed stars in NGC 2419. 
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